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I .  Financial  Report . 

The  University  of  Wisconsin  Modelling,  Information  Processing  and  Control  (MIPAC) 
Facility  began  operation  in  October,  1983,  with  initial  acquisitions  of  equipment  provided 
for  under  the  OoO  University  instrumentation  Program.  The  equipment  was  funded  by  two 
separate  grants,  one  from  the  Air  Force  Office  of  Scientific  Research  (Grant  No.  AFOSR-83- 
0281,  $110,000)  and  one  from  the  Army  Research  Office  (Grant  No.  DAAG29-83-G-0056, 

$28,000) .  The  following  items  were  acquired  from  Hewlett-Packard,  Inc.  under  Grant  No. 

AF0SR-83-0281  in  accordance  with  University  of  Wisconsin  policies  and  procedures. 


Cost 


1.  Purchase  Order  No.  UAD  078Y002  Dated:  9/29/83  $101,706.01 

This  purchase  order  covered  the  initial  acquisition  of  our 
Hewlett  Packard  5451C  Fourier  Analyzer  and  related  peripheral 
equipment.  Specifically,  the  items  purchased  from  Hewlett- 
Packard,  Inc.  are  as  follows: 

5451C  Fourier  Analyzer,  Options  600  &  S23 
10833B  HPIB  Cable 

10697D  Field  Add-on:  S83  Var  Resis  Scope 

620  DAC  Module 
350  Vibration  Con 
451  Less  5477A 

After  some  months  of  initial  operation  it  became  clear  that 
more  flexible  use  of  the  Hewlett-Packard  5451C  Fourier 
Analyzer  required  acquisitions  and  installation  of  a  different 
operating  system  for  the  attached  HP  1000  computer  and  a  tape 
cassette  unit  in  the  Fourier  Analyzer's  HP  2648A  graphics 
terminal.  The  complete  acquisition  from  Hewlett-Packard, 

Inc.  which  covered  these,  and  related,  items  is  summarized  below: 

2.  Purchase  Order  No.  UAD  924C055  Dated;  6/30/84  6395.50 

12747H  128  KByte  Memory  Board 
92068A  RTE  IVB  System 
92068A  Opt.  031,  7906  Disc  Media 
12539C  Time  Base  Generator 

12966A  Buffered  Asynchronous  Interface,  Option  001 
13236B  Cartridge  Tape  Refit  for  HP2648A  Graphics  Terminal 
13261A  2648A  ROMS  (Software  for  Cartridge  Tape  Unit)  ITR  252 

The  final  acquisitions  under  the  subject  grant  were  three  disk 
cartridges  and  four  boxes  of  tape  cassetes  for  use  with  the 
HP  5451C's  hard  disk  unit  and  tape  cassette  unit,  respectively. 

These  are  listed  below. 

3.  Purchase  Order  No.  UAE  924C545  Dated:  9/27/84  980.00 

5  Cartons  HP  98200  A  Cartridge  Tapes 

4  HP  12940  A  Disc  Cartridges  _ 


TOTAL 


$109,081.51 


MODELLING  AND  PARAMETER  ESTIMATION 


II.  Scientific  Report 


FOR  DISTRIBUTED  VIBRATORY  SYSTEMS;* 


A  Report  on  Research  Projects  at  the  MIPAC  Facility 
of  the  University  of  Wisconsin/  Madison/  to  June  30/  1984 


by  D.  L.  Russell** 


ABSTRACT 

Citing  data  collected  and  analyzed  with  the  University  of  Wisconsin  MIPAC 
(Modelling/  Information  Processing  and  Control)  Facility  electronic  measurement  and 
analysis  equipment,  we  make  some  general  comments  concerning  mathematical  models 
which  appear  to  be  appropriate  for  modelling  certain  vibratory  systems  of 
distributed  parameter  type.  Aspects  discussed  include:  location  of  vibrational 
spectra,  damping  rates,  and  spectral  displacement  due  to  ma?  t.  density  and/or 
elasticity  variations.  Particular  emphasis  is  placed  on  some  properties  of 
segmented  beams.  The  article  ends  with  a  preliminary  mathematical  discussion  of  the 
feasibility  of  parameter  identification,  from  vibrational  spectrum  data  alone,  in 
the  wave  and  Euler-Bernoulli  beam  equations. 


♦Supported  through  the  DoD  University  Research  Instrumentation  Program  (FY  83):  Air 
Force  Office  of  Scientific  Research  under  Grant  No.  AFOSR-83-0281  and  Army  Research 
Office  under  Grant  No.  DAAG29-83-G-0056. 


♦♦Mathematics  Research  Center,  University  of  Wisconsin,  Madison,  WI  53706,  U.S.A. 


MODELLING  AND  PARAMETER  ESTIMATION 


FOR  DISTRIBUTED  VIBRATORY  SYSTEMS 

1.  Introduction. 

Beginning  in  October,  1983,  we  have  been  developing  the  University  of 
Wisconsin  Modelling,  Information  Processing  and  Control  (MIPAC)  Facility 
with  the  cooperation  of  the  Department  of  Mathematics,  the  Department  of 
Electrical  and  Computer  Engineering  and  the  Mathematics  Research  Center.  Our 
objectives  in  the  operation  of  this  Facility  are  as  follows: 

1.  Better  understanding  of  physical  processes: 

(a)  From  the  dynamical  point  of  view; 

(b)  From  the  point  of  view  of  measurement  capabilities  and  limitations; 

(c)  From  the  point  of  view  of  our  ability  to  model  the  process  mathema¬ 

tically. 

2.  Development  of  modelling  techniques: 

(a)  Qualitative  model  identification; 

(b)  Model  calibration  and  parameter  estimation. 

3.  Development  of  information  processing  techniques: 

(a)  Instrument  capabilities  and  limitations; 

(b)  Noise  limitations; 

(c)  Discretization  problems; 

(d)  Mathematical  information  processing. 

4.  Development  of  real  time  control  techniques: 

(a)  Actuator  capabilities  and  limitations; 

(b)  Real  time  computational  capabilities  and  limitations,  including 

discretization  problems; 

(c)  Mathematical  control  law  development. 

Initial  grants  from  the  Air  Force  Office  of  Scientific  Research  and  the  Army 


Research  Office,  with  additional  assistance  from  the  Wisconsin  Alumni  Re¬ 
search  Foundation,  have  enabled  us  to  acquire  and  begin  operation  of  a  two 
unit  Facility,  incorporating  an  Analysis  Unit  and  a  Model  Development  Unit. 

The  electronic  equipment  acquired  includes  a  Hewlett-Packard  5451C  Fourier 
Analyzer,  a  Bruel  §  Kjaer  Vibration  Exciter  and  supporting  peripheral  units. 

The  University  is  providing  us  with  approximately  1500  square  feet  of  newly 
remodelled  laboratory  space  which  we  expect  to  overcrowd  very  quickly.  At 
present  two  faculty  members  and  two  research  assistants  are  working  with  the 
facility  and  we  expect  to  add  an  additional  senior  staff  member  in  the  autumn 
of  1984. 

The  paragraphs  and  figures  which  follow  do  not  constitute  a  research 
paper  as  such  but,  rather,  a  description  of  the  work  which  we  have  been  carry¬ 
ing  out  in  the  MIPAC  Facility,  indications  for  future  modelling  work,  both  of 
a  mathematical  character  and  of  the  laboratory  sort,  and  sketches  of  some  of 
the  mathematical  questions  likely  to  enter  into  our  work.  Not  everything  is 
reported  here;  we  have  not,  for  example,  discussed  our  extensive  laboratory 
work  with  various  types  of  nonlinear  oscillators . 

One  point  which  we  wish  to  make  very  strongly  is  the  way  in  which  physi¬ 
cal  data  can  prompt  and  direct  mathematical  investigation  in  unexpectedly  signi¬ 
ficant  ways.  We  will  illustrate  this  in  Section  3  by  exhibiting  data  obtained 
from  segmented  beams  and  discussing,  briefly,  the  implications  as  regards  the 
mathematical  work  required  to  provide  an  anlytical  framework  to  correspond  to 
and  explain  what  we  observe. 

In  Section  2  of  this  article  we  discuss  data  obtained  from  vibrations  of 
uniform  beams.  This  provides  a  necessary  background  for  Section  3  and  also 
allows  us  to  discuss  important  laboratory  evidence  in  regard  to  the  relation¬ 
ship  between  the  frequency  of  a  mode  of  vibration  and  the  rate  at  which  that 


mode  of  vibration  is  attenuated  through  internal  energy  decay. 

.  The  last  section  of  this  article,  Section  4,  is  an  analysis  of  the  feasi¬ 
bility  of  first  order  mass  density  estimation/identification  on  the  sole  basis 
of  observed  spectral  data.  We  will  see  how  a  very  practical  engineering 
problem  leads  to  some  quite  interesting  and  intricate  mathematical  questions 
connected  with  completeness  and  linear  independence  of  certain  sets  of  functions 
derived  from  the  eigenfunctions  of  a  vibrating  system. 

We  will  freely  use  terminology  common  in  the  literature  of  discrete 
Fourier  methods.  The  reader  is  referred  to  I  5  ]  for  definition  of  these 
terms  and  development  of  the  related  theory. 


-vr'r 


2.  Frequency  Distribution  and  Damping  for  Uniform  Thin  Beams. 

There  are  a  number  of  ways  in  t;hich  the  natural  frequencies  of  vibration 
of  an  elastic  body  may  be  determined.  The  one  which  we,  along  with  others, 
have  used  in  connection  with  various  types  of  beams  is  extremely  simple  but 
quite  effective.  The  beam  is  tapped,  once,  and  quite  smartly,  with  an  appropri¬ 
ate  hammer  or  similar  instrument.  The  resulting  vibrations  of  the  beam  are 
transformed  into  electric  signals  by  an  accelerometer  attached  to  the  beam. 

The  data  acquisition  functions  of  the  Fourier  analyzer  are  triggered  and  the 
accelerometer  signals  are  read  in  and  recorded.  A  record  such  as  is  shown  in 
Figure  1  is  obtained.  The  corresponding  discrete  Fourier  transform  and  the 
associated  power  spectrum  are  computed.  This  procedure  is  repeated  a  number 
of  times,  typically  25  to  50,  and  an  averaged  power  spectrum  is  obtained, 
whose  logarithm  is  then  computed  and  displayed,  as  in  Figure  2.  In  the  process 
the  beam  is  struck  at  random  points  along  its  length  to  ensure  excitation  of  all 
modes . 

The  log  averaged  power  spectrum  (LAPS)  for  a  thin  beam  (approximately  2  mm. 

thick,  25  ram.  wide  and  1  meter  long)  of  mild  steel,  clamped  at  one  end  and  free 

at  the  other,  obtained  as  indicated  above,  is  shown  in  Figure  2  for  the  frequency 

* 

range  0  -  2500  Hz.  Depending  on  where  one  stops  counting,  11  or  12  distinct  modes 
of  vibration,  whose  associated  natural  frequencies  are  indicated  by  the  sharp 
peaks  in  the  LAPS,  appear.  Beyond  about  1200  Hz.  the  peaks  become  indisting¬ 
uishable  from  the  ambient  noise  and  are  no  longer  significant.  The  measurement 
used  in  this  case  was  the  acceleration  of  the  beam  at  the  free  boundary. 

In  Figure  3  we  show  the  same  spectrum,  now  restricted  to  the  range  0  -  1000  Hz., 
in  ordex>  to  compare  it  with  similar  LAPS  obtained  from  other  types  of  beams.  Fig¬ 
ure  4  exhibits  the  spectrum  from  a  thin  aluminum  beam.  This  beam  is  highly  flex¬ 
ible  and,  as  a  result,  the  number  of  spectral  peaks,  corresponding  to  natural 
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frequencies,  occurring  in  the  range  0  -  500  Hz.  here  is  comparable  to  the  num¬ 
ber  occurring  in  twice  that  range  in  Figure  3.  In  Figure  5  we  see  the  LAPS 
obtained  from  a  thin  hardwood  beam.  In  this  case  a  much  smaller  number  of  nat¬ 
ural  modes  of  vibration  is  observable,  perhaps  because  higher  modes  are  over¬ 
damped  in  this  case.  In  Figure  6  we  see  the  spectrum  obtained  from  a  segmented 
wooden  beam,  actually  a  good  quality  carpenter's  rule  obtained  at  a  local  hard¬ 
ware  store.  The  most  significant  feature  occurring  in  this  case,  the  gap  bet¬ 
ween  approximately  140  Hz.  and  280  Hz.,  will  be  discussed  in  Section  3  below. 

The  decay  of  amplitude  with  increasing  frequency  evident  in  Figures  4 
through  6,  in  contrast  to  Figure  3,"  is  not  particularly  significant.  Figure  3 
was  obtained  using  acceleration  data  while  Figures  4  through  6  were  obtained 
from  velocity  data. 

Since  we  are,  eventually,  interested  in  control  applications,  we  will  place 
particular  emphasis  on  the  properties  of  systems  which  are  important  in  this  con¬ 
nection.  A  very  old  argument,  re-disputed  by  every  new  generation  of  control 
theorists,  concerns  the  necessity  and/or  advisability  of  distributed  parameter 
representation  of  spatially  extended  vibrating  continua  such  as  beams,  plates, 
shells,  etc.  We  are  hardly  prepared  to  pronounce  judgment  on  this  unresolvable 
question  but  we  are  in  a  position  to  make  a  few  relevant  remarks. 

First  of  all,  it  is  evident  that  any  reasonable  answer  to  the  question  must 
take  into  account  the  mode  of  anticipated  excitation  and  the  mode  of  output  evalu¬ 
ation.  Here  we  restrict  attention  to  the  latter.  Figure  3  corresponds  to  accele¬ 
ration  measurements.  Here  the  higher  frequencies  do  not  fall  off  very  markedly, 
in  contrast  to  Figures  4  through  6,  taken  with  velocity  measurements  where,  if 
we  take  into  account  that  this  is  a  logarithmic  plot,  nothing  beyond  the  first 
two  or  three  modes  is  of  much  significance. 

If  we  assume  the  worst  case  (best  case  for  the  proponents  of  distributed 


parameter  theory)  corresponding  to  Figure  3,  we  see  that  adequate  modal  descrip¬ 
tion  of  the  system  might  require  use  of  a  linear  system  whose  dimension  is  on 
the  order  of  20.  A  distributed  parameter  approach  necessarily  presupposes 
the  use  of  some  approximation  procedure,  e.g.,  one  based  on  spline  functions, 
in  the  final  application.  Assuming  approximation  via  cubic  splines,  a  system 
dimension  of  20  allows  subdivision  of  the  spatial  domain  into  about  8  sub inter¬ 
vals.  Comparison  should  be  made  between  modelling  and  control  results  obtained 
with  this  number  of  subintervals  and  corresponding  results  obtained  using  a 
modal  representation.  The  frequencies  present  in  the  model  will  not  be  identi¬ 
cal  in  the  two  approaches  and,  perhaps  more  importantly,  the  control  and/or 
disturbance  input  coefficients  (£  61)  obtained  via  the  two  approaches  will 
differ.  We  feel  that  it  will  be  interesting  to  see  how  the  performances  compare 
when  the  two  contenders  are  matched  to  each  other  in  this  way. 

Without  pre- judging  the  outcome  of  the  contest  we  have  just  proposed,  let  us 

proceed  under  the  assumption  that  a  distributed  parameter  model  is  to  be  used. 

The  simplest  thin  beam  model  based  on  partial  differential  equations  is  provided 

by  the  Euler-Bernoulli  equation 

Pw  +  El  w  ■  0, 
tt  xxxx 

where  is  the  mass  per  unit  length  and  El  is  the  bending  modulus.  In  the  un¬ 
forced,  clamped  case,  which  corresponds  to  our  data,  the  appropriate  boundary  con¬ 
ditions  are 

w(0,t)  -  0,  w  (0,  )  *  0, 
x 

W  (l,t)  a  0,  w  (L,t)  «  0, 

XX  XXX 

where  L  is  the  length  of  the  beam.  With 

iwt 

w(x,t)  a  e  z(x) 


we  find  that 


% 

and  with  oi  «  (  p/EI)  we  must  have 

'  h  H  h  h 

z(x)  ■  c  sin(o<W  x)+c  cos(e<td  x)+c  sinh(tf&?  x)+c  cosh(<rt«o  x) . 
12  3  4 

The  boundary  conditions  at  x  »  0  imply 

c  ♦  c  *  0,  c  +  c  *  0, 

2  4  13 

so  that,  with  c  *  -c  »  c,  c  «  -c  »  d, 

13  2  4 

Jj 

z(x)  »  c(sin(tf&>  x)  -  sinhCtfto  x))  +  d(cos(<*<*>  x)  -  cosh(  x)) . 
Taking  the  second  and  third  derivatives  at  x  *  L  we  obtain  the  equations 

c(-sin(«^  L)-sinh(«<w  iOJ+dC-cosC*^  L)-cosh( l))  *  0, 

h  h  h  h 

c(-cos(rt«  L)-cosh(<*<*>  L))+d(  sin(*&>  L)-sinh(*C4>  L))  *  0. 

The  obvious  determinant  condition  then  gives 

jj  jj  jj 

(-cos(o(k>  L)-cosh(o(  u>  L))C~cos(^^'  Lj-coshC  l)) 

h  h  H  h 

-(sin(«<<*>  L)+sinh(«*  *<-’  L))(-sin(<rfc^  L)+sinh(^  &  L))  *  0 
or,  taking  familiar  identities  into  account, 

H  h 

1  ♦  cos(c(cO  L)cosh(*w  L)  »  0. 

Examination  of  the  data  from  which  Figure  2  was  obtained  reveals  the  fol- 

% 

lowing  approximate  values  for  the  frequencies  CO  and  their  square  roots  (Co  ) 

k  k 


*k 

tco  ) 
k 

2.97 

1.72 

19.54 

4.42 

53.72 

7.33 

107.50 

10.37 

178. 

13.34 

266. 

16.31 

374. 

19.34 

495.5 

22.26 

635. 

25.20 

788.5 

28.08 

942. 

30.69 

Since,  asymptotically,  ©<(CJ  )  L>  (k  -  *s)  ,  we  obtain  the  approximate 

k 

value  which  must  be  assigned  to  XL  (from  averaging  the  last  five  values  in 


the  table) 


XL  *  (1/S)  | 
«  1.062. 


/ 10. 

^30. 


S_  .  9.5  +  8.5  •  7.5  6.5  \ 

69  28.08  25.20  22.26  19.34/ 


using  this  value  to  predict  CO  via  solution  of  the  equation  provided  by  the 

1 

determinant  condition  we  obtain  the  value 

( xi  *  1.76 
1 

which  may  be  compared  with  the  observed  value  1.72.  As  the  observed  value  may 
be  significantly  lowered  by  movement  of  the  supports  and  the  mass  of  the  acceler¬ 
ometer,  we  may  regard  this  as  fairly  good  agreement.  Thus  the  Euler-Bemoulli 
model  appears  to  be  reasonably  accurate  in  predicting  the  observed  frequencies 
of  the  beam  in  this  case. 

This  example  does  point  out,  however,  one  aspect  of  the  distributed  parameter 
approach  which  must  be  taken  into  account.  This  approach  presupposes  a  certain 
model  or  class  of  models.  If  exact  spectral  matching  is  important  it  may  be  a 
serious  disadvantage  that  no  member  of  the  model  class  is  actually  able  to  achieve 


^  '  ■  >l  J. • turn  u.-  -_B  LHlWUi'iLIJmT*.!! 


the  observed  spectral  pattern.  This  is  true  for  the  observed  beam  frequencies 
shown  above  if  the  Euler- Bernoulli,  model  is  assumed  at  the  cutset:  no  spatially 
constant  ^  ,  El  can  achieve  exactly  the  spectral  pattern  shown,  even  though 
the  beam  is,  to  every  outward  appearance,  completely  uniform.  We  will  have  more 
to  say  about  variable  coefficients  in  Section  4  of  this  article. 

It  seems  likely  that  a  distributed  parameter  model,  such  as  the  Euler-Ber- 
noulli  model  or  some  finite  element  approximation  thereof,  will  have  significant 
advantages  from  the  point  of  view  of  on-line  identification  because  of  the  sna  Her 
number  of  parameters  necessary  in  order  to  describe  the  model,  as  compared,  for 
example,  with  a  set  of  equations 


x  «-Xx  +bu  (u  *  the  applied  control) 
k  k  k  k 


derived,  from  a  modal  description.  .  While  the  latter,  because,  of  the  large,  number 
of  parameters  involved,  permits  more  or  less  exact  frequency  matching,  the  very 
multiplicity  of  parameters  present  makes  for  unstable,  ill-conditioned  identi¬ 
fication.  The  more  modest  number  of  parameters  involved  in  a  model  with  an  a 
priori  assumed  structure  makes  the  identification  task  a  much  moTe  tractab!  i 
problem. 

The  Euler-Bemoulli  equation  in  its  commonly  written  form 

jew  +  EIw  ■  0 
'  tt  xxxx 

corresponds  to  conservation  of  the  energy  form 


*  fL  2  2 

C  9  h]  C(0(w  )  +  EI(w  )  )dx, 

‘'O  V  .  t  xx 


where  w  is  lateral  deflection  and  L  is  the  length  of  the  beam.  The  constants 
^  and  El  are  mass  per  unit  length  and  bending  modulus,  respectively.  It  is  well 
known  that  energy  losses  in  fact  occur  and  that  these  losses  are  frequency  depend¬ 
ent.  This  is  amply  documented  by  inspection  of  Figure  1  wherein  it  is  clear  that 
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the  high  frequency  components  of  the  motion  are  attentuated  much  more  rapidly 
than  are  the  low  frequency  components. 

A  quantitative  appreciation  of  the  dependence  of  damping  rate  on  frequency 
can  be  obtained  by  taking  a  number  of  consecutive  time  records  after  each  exci¬ 
tation  of  the  beam  and  computing  the  corresponding  power  spectra  separately. 

The  LAPS  may  then  be  overlayed  and  the  vertical  positions  of  the  peaks  at  the 
same  frequency  compared  in  order  to  estimate  the  damping  rate.  In  Figure  7, 
which  shows  the  fr^.  .  r.cy  r;.r.ge  0  -  500  Hz.  for  the  same  beam  as  was  used  in 
Figures  2  and  3,  -our  such  consecutive  spectra  are  overlayed,  the  peaks  being  ; 
darkened  for  easy  comparison.  The  vertical  gap  between  successive  peaks  is 
proportional  to  the  damping  coefficient  for  the  mode  with  the  indicated  fre¬ 
quency.  Our  experiments  indicate  a  rough  linear  proportionality  between  the 
damping  rate  and  the  frequency  in  modes  1,2,4  and  5,. becoming  tore  nearly  con¬ 
stant  in  modes  5  through  7.  The  anomalous  behavior  of  mode  3  is  thought  to  be 
due  to  energy  transfer  to  a  torsional  vibration  mode  whose  frequency  is  close 
to  that  of  the  third  lateral  mode  shown.  The  linear  dependence  of  damping  rate 
on  frequency  is  more  impressive  in  Figure  8,  which  shows  tuo  successive  overlayed 
spectra  for  a  different  beam. 

A  linear  damping  /  frequency  ratio  corresponds  to  what  is  called  struc¬ 
tural  damping  in  the  engineering  literature  I  4  J .  As  shown  in  [  3  1 ,  it  is 
consistent  with  the  model 

fw  +#Aw  +  Aw  »  0 

tt  t 

h 

where  A  is  the  operator  EIw  defined  on  an  appropriate  domain  and  A  is  the 

xxxx 

unique  positive  definite  square  root  of  A.  It  is  also  consistent  with  computed 

exponential  solutions  of  the  partial  differential  equation 

*>w  -  ^w  ♦  EIw  *  0 

V  tt  txx  xxxx 
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with  the  boundary  conditions 

w(o,t)  *  0,  w  (0,t)  *  0, 
x 

EIw  CL,t)-yw  (L,t)  -  0,  w  (L,t)  »  0. 

XXX  tx  XX 

The  semigroup  theory  related  to  this  equation  may  be  based  on  the  energy  dissi¬ 
pation  law 

(d/dt)  S  (w(\t),w  (*,t))  *  -Jf/  (w  (x,t))  dx£0 
t  v  0  tx 

but  it  remains  to  be  shown  that  the  generator  is  a  spectral  operator,  which  ist 
necessary  if  we  are  to  be  able  to  draw  conclusions  from  the  spectrum  of  the 
operator,  and  that  the  semigroup  is  a  holomorphic  semigroup.  It  is  clear -that 
this  system  deserves  substantially  more  attention  in  the  future. 
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3.  Vibration  Spectra  of  Segmented  Beams. 

Most  large  structures,  including  the  rather  exotic  examples  that  have  come 
to  be  known  as  "large  space  structures",  are  not  single  units  but  are,  rather, 
assemblages  of  many  small  units  joined  together  in  various  ways.  Thus,  in  many 
respects,  simple  systems  such  as  the  wave  equation,  Euler-Bernoulli  beam  equa¬ 
tion,  thin  plate  quation,  etc  ,  have  unrealistically  simple  structure  as  com¬ 
pared  with  real  systems.  On  the  other  hand  it  is  clear  that  one  cannot  develop 
a  theory  of  distributed  parameter  systems  on  the  basis  of  structures  of  arbi¬ 
trary  complexity;  there  has  to  be  some  unifying  pattern  and,  inevitably,  this 
involves  a  degree  of  simplification.  One  is  led  to  seek  out,  as  a  consequence, 
a  class  of  systems  simple  enough  to  admit  some  hope  of  analytical  treatment  and 
yet  possessing  enough  intricacy  to  carry  one  beyond  the  very  simplest  mod  s. 

It  seems  to  us  that  the  models  representing  the  vibration  of  segmented  beams  are 
ideally  situated  on  this  middle  ground  between  simplicity  and  complexity. 

For  the  purposes  of  this  article  we  will  consider  a  beam  occupying  a  spatial 
interval  which,  without  loss  of  generality,  we  may  take  to  be  [0,Tf  I ,  and  simply 
supported  at  x  *  0  and  x  =  Tf  .  By  "simply  supported"  we  mean  that  the  boundary 
conditions 

w(0,t)  =  0,  w  (0,t)  *  0,  w(tf,t)  *  0,  w  (tf, t)  *  0 

XX  XX 

apply. 

We  will  suppose  that  the  interval  [0,1^  ]  is  divided  into  N  segments  of  eq¬ 
ual  length  by  means  of  the  points 

x  *  0,  x  =  (klV)/N,  k  =  1,2,...,N. 
o  k 

In  each  of  the  intervals  I  «  lx  ,x  I,  k  ■  1,2,...,N,  we  assume  that  the  Euler- 

k  k-1  k 

Bernoulli  beam  equation  (here  simplified  with  unit  coefficients) 
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is  satisfied  (whether  in  a  classical  or  distribution  sense  need  not  concern 


us  here)  while,  at  the  junction  points,  x  ,  we  impose  conditions  which  reflect 

k 

our  assumptions  about  the  nature  of  the  joints  there.  Here  we  will  distinguish 
only  two  different  types,  shown  schematically  in  Figure  13: 


Flexible  joints:  in  this  case  the  energy  expression  is 

rtf  2  2  N  2 

£  *  h  ((w  (x,t)  +  (w  (x,t))  )dx  +  %iS>Cw  (x  +,t)-w  (x  -,t))  , 

v  0  t  xx  k»l  x  k  x  k 

^denoting  the  stiffness  coefficient  for  the  joints.  Conservation  of  energy 

requirements  then  lead  to  the  partial  differential  equation  indicated  above 

and  the  junction  conditions 

v  (x  :,t)  =  w  (x  -,t)  *  ff'Cw  (x  -  w  (x  -,t)),  - 
xxk  xxk  xk  xk 

w  (x  +,t)  =  w  (x  -,t), 
xxx  k  xxx  k 

w(x  +,t)  =  w(x  -,t),  k  *  1,2, . . .,N-1. 
k  k 

Here  f(x+)  means  lim  f(x+y),  etc. 
ylO 

Massive  stiff  joints:  here  we  suppose  that  no  macroscopic  degree  of  bending 

occurs  at  the  joints,  which  are  massively  reinforced  by  structures  having  very 

little  spatial  extent  but  appreciable  mass.  In  this  case,  if  the  reinforcing 

masses  have  mass  m,  the  energy  expression  becomes 

»  ftf  2  2  N  2 

O  *  h  ((w  (x,t))  +  (w  (x,t))  )dx  +  Jsm^(w  (x  ,t))  . 
vO  t  xx  k=l  t  k 

Energy  conservation  assumptions  then  lead  to  the  same  partial  differential  eq¬ 
uations  as  before  and  the  junction  conditions 


w(x  +,t)  =  w(x  -,t), 
k  k 

w  (x  +,t)  =  w  (x  -,t),  w  (x  +,t)  ■  w  (x  -,t), 
xk  xk  xxk  xx  k 

mw  (x  ,t)  =  w  (x  -,t)  -  w  (x  +,t),  k  =  1,2,...,N-1. 
tt  k  xxx  k  xxx  k 

Before  discussing  the  theoretical  aspects  of  these  systems,  let  us  consider 
some  of  the  experimental  results .  These  were  obtained  with  a  segmented  wooden 
beam,  consisting  of  short,  thin  hardwood  segments,  fastened  to  each  other,  end 
to  end  (a  folding  carpenter's  rule  was  used,  in  fact).  The  resulting  segmented 
beam  was  clamped  at  one  end  and  left  free  at  the  other;  thus,  as  far  as  the  boun¬ 
dary  conditions  are  concerned,  the  experimental  apparatus  does  not  exactly  match 
the  theoretical  model  described  above.  However,  we  expect  the  effect  of  different 
boundary  conditions  to  be  small  and  quantitative  rather  than  qualitative.  The 
experimental  beam  has  slightly  flexible  joints,  including  the  joint  at  the  point 
where  it  is  clamped. 

In  different  experiments,  varying  numbers  of  segments  of  the  beam  were  ex¬ 
tended  out  from  the  clamping  apparatus.  In  each  case  the  beam  was  excited  by 
striking  and  the  LAPS  was  obtained  using  the  Fourier  analyzer.  The  resulting 
power  spectra  for  beams  composed  of  5,  7  and  8  projecting  segments  are  shown  in 
Figures  9  through  11,  respectively.  It  will  be  observed  in  each  case  that  S,  7 
and  8  vibrational  frequencies  appear  at  the  low  end  of  the  spectrum  and  then  there 
is  a  gap  until  the  next  frequency  appears,  as  indicated  by  the  peak  just  to  the 
left  of  the  300  Hz.  point.  Another  group  of  frequencies  occurs,  not  too  dis¬ 
tinctly,  from  300  Hz.  to  400  Hz.,  approximately. 

In  Figure  12  the  preceding  three  figures  are  superimposed.  What  is  striking 
is  the  fact  that  the  first  group  of  5,7  and  8  frequencies,  depending  on  the  num¬ 
ber  of  beam  segments,  fills  approximately  the  same  frequency  range  in  each  case. 
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Figure  9 


Figure  10 


Figure  11 


while  the  next  frequency  appears  at  roughly  the  same  location  in  each  case. 

How  is  this  pattern  to  be  accounted  for?  We  will  see  that  examination  of  the 
theorectical  model  provides  an  explanation.  We  will  do  this  only  in  the  context 
of  the  flexible  joint  model  but  we  will  make  some  remarks  relative  to  the  massive 
stiff  jointed  model  as  well. 

Let  us  note  that  in  the  simply  supported,  N  segment  case,  if  the  segments 
were  infinitely  stiff  by  comparison  with  the  joints,  we  would  expect  to  find  a 
certain  number  of  modes  of  vibration  corresponding  to  bending  of  the  structure 
at  the  joints,  the  segments  remaining  straight..  It  is  easy  to  see  that  the  re¬ 
sulting  linear  dynamical  system  has  dimension  2(N-1);  there  will  exist  N-l  os¬ 
cillatory  modes  corresponding  to  bending  at  the  joints.  Ac  the  stiffness  of  the 
beam  segments  decreases  relative  to  that  of  the  joints,  these  vibrational  fre¬ 
quencies  will  decrease  and  the  mode  shapes  will  change  somewhat  but  the  character 
of  these  first  N-l  modes  will  remain  much  the  same  as  long  as  the  beam  segments 
remain  quite  stiff  relative  to  the  joints. 

Now  let  us  consider  the  N-th  mode.  In  the  unsegmented,  uniform  case  it  is 
easy  to  see  that  the  modal  function,  i.e.,  eigenfunction,  in  this  case  is 

0  (x)  *  sin(Nx) . 

N 

However,  because  this  function  has  zero  values  and  second  derivatives  at  the 

points  x  ,  k  *  1,2,..., N-l,  all  of  the  boundary  conditions  there  imposed  on  the 
k 

segmented  beam  are  satisfied  by  0  (x)  and  we  conclude  that  this  is  the  N-th 

N- 

eigenfunction  in  the  segmented  case  as  well  -  no  matter  what  the  value  of  N  may 
be.  This  reasoning  will  apply  to  the  2N-th,  3N-th,  etc.,  modes  as  well.  One  can 
also  see  that  the  frequency  of  vibration  of  the  mN-th  mode  agrees  with  the  m-th 
modal  frequency  of  a  beam  of  length  Xf /  N .  If  we  pose  the  same  problem  on  an 
interval  of  length  NIT'  ,  which  corresponds  to  the  experimental  reality  where 
the  length  of  the  beam  depends  on  the  number  of  projecting  segments,  we  then  see 
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that  the  frequency  of  vibration  of  the  mN-th  mode  is  invariant  -  it  always  agrees 

with  that  of  the  m-th  mode  of  vibration  of  a  simply  supported  beam  of  length  It  . 

This  observation  corresponds  to  the  near  invariance  of  the  (N+l)-st  mode  in  our 

experiments  (the  index  shifts  from  N  to  N+l,  however,  because  the  experimental 

beam  is  free  at  one  end  and  a  flexible  joint  exists  at  the  clamped  end,  thus 

introducing  an  additional  mode  of  vibration  for  the  infinitely  stiff  case) . 

In  the  uniform,  unsegmented  case,  the  first  N-l  frequencies  of  the  simply 

supported  beam  would  be  distributed  throughout  the  interval  between  0  and  the 

frequency  of  the  N-th  mode,  the  k-th  modal  frequency  being  roughly  proportional 
2  ' 

to  k  .  Because  of  the  additional  flexibility  permitted  by  bending  at  the  joints, 

however,  those  frequencies  are  all  significantly  lowered  in  the  segmented  case 

(recall  that  the  k-th  eigenvalue  of  the  operator  EIw  minimizes  the  so-called 

xxxx 

Rayleigh  quotient 

jrv  2  / r*  2 

(w  )  dx  /  /  w  dx 
0  xx  /  J  0 

subject  to  the  imposed  boundary  conditions  and  orthogonality  to  the  first  k-1 
eigenfunctions;  additional  flexibility  results  in  a  numerator  which  is  somewhat 
Sjnalier  relative  to  the  denominator.  Since  the  N-th  modal  frequency  is  unchanged 
as  we  pass  from  the  uniform  to  the  segmented  case,  a  "gap”  opens  up  between  the 
N-th  mode  and  the  (N-l)-st  mode.  (This  corresponds  to  the  gap  between  the  N-th 
and  (N+l) -st  modes  in  the  experimental  situation  which  we  have  described.) 

Similar  gaps  may  be  expected  to  occur  between  the  mN-th  and  (mN-l)-st  frequen¬ 
cies  with  the  (m-l)N-th  through  (mN-l)-st  frequencies  being,  relatively  speaking, 
clustered  together. 

Examination  of  the  second  type  of  system,  involving  massive,  stiff  joints, 
may  be  expected  to  uncover  a  similar  pattern  of  clusters  and  gaps;  again  the  mN-th . 
modes  will  not  change  from  the  uniform  case,  the  massive  joints  remaining  motion¬ 
less  at  the  nodes,  while  the  intermediate  modes  will  exhibit  depressed  frequeifies 


due  to  the  additional  inertial  effects  of  those  masses  (which  make  the  denominator 
of  the  Rayleigh  quotient  larger) . 

The  results  sketched  here  may  be  expected  to  have  some  importance  in  control 
applications .  The  spectral  gaps  which  we  have  described  can  be  used  to  advantage 
to  permit  control  and/or  observation  of  the  clustered  modes  in  the  first  group 
without  excessive  interference  from  the  modes  in  the  second,  or  higher,  groups 
which  will  be  separated  from  the  first  group  by  the  spectral  gap  whose  origins 
we  have  discussed  here.  It  is  clear  that  the  subject  deserves  more  careful  and 
precise  mathematical  treatment,  incorporating  estimates  and  asymptotic  expressions 
for  the  eigenvalues  and  eigenvectors  of  these  segmented  beam  systems.  In  this 
connection  we  are  happy  to  note  that  some  contributions  in  this  direction  have 
recently  been  made  by  Chen  121. 
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4.-  Parameter  Estimates  from  Frequency  Data. 


The  subject  matter  of  this  section  is  best  illustrated  by  reference  to 
Figure  14.  The  thin  steel  beam  from  which  the  log-averaged  power  spectra  of 
Figures  2  and  3  were  obtained  was  first  excited  in  its  original  configuration 
and  then,  after  its  mass  had  been  modified  by  the  addition  of  several  small 
permanent  magnets  near  its  center,  as  illustrated  in  Figure  15.  The  two 
power  spectra  were  superimposed  to  produce  Figure  14.  In  this  figure  the  super¬ 
imposed  power  spectra  produce  double  peaks,  i.e.,  peaks  lying  near  each  other. 

In  each  case  the  peak  to  the  left  corresponds  to  the  beam  with  additional  mass, 
as  would  be  expected.  A  very  natural  question  then  arises:  is  it  possible  to 
reconstruct  the  modified  mass  distrioution  from  the  observed  perturbations  of 
the  vibration  spectrum?  Our  purpose  in  this  section  of  the  paper  is  to  discuss 
a  first  order  approach  to  this  problem  from  a  mathematical  point  of  view.  Some 
theoretical  studies  yielding  more  or  less  global  uniqueness  results  for  mass/ 
elasticity  distributions  versus  frequency  spectra  obtained  with  various  boundary 
conditions  have  been  presented  by  McLaughlin  1 10  1 ,  £  11  ] . 

Let  us  consider  a  second  order  linear  system  in  a  Hilbert  space  X: 


y  ♦  Ay  *  0,  ye«#(A)  C  x, 

(densely) 


where  A  is  an  unbounded  self-adjoint  operator  on  X  with  eigenvalues  X  ,  k  »  1, 

k 

2,3,...  of  finite  multiplicity  and,  without  loss  of  generality,  nondecreasing. 


with 


lim  X  *  oo, 
k*oo  k 


the  corresponding  normalized  eigenvectors  0  ,  k  *  1,2,3,...,  being  selected  so 

k 

as  to  form  an  orthonormal  basis  for  X.  We  let 


-/T .  X 


CJ  “VA  >  £ 

k  k  k 
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Figure  14 


(a)  Clamped  /  Free  Beam  in  Original  Configuration, 


x  ■  0  x  «  L 


(b)  Clamped  /  Free  Beam  with  Added  Permanent -Magnet  Masses. 
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As  is  well  known,  the  above  second  order  equation  is  associated  with  the 


first  order  system 


*(:)•(-:  !)(:) 


(x  «  y) 


and  the  operator 


A  "U  l) 


generates  a  strongly  continuous  group  of  bounded  operators  on  Y  x  X,  where  Y  is 

the  closure  of  o^(A)C  X  with  respect  to  ((  ||  ,  the  norm  associated  with  the  inner 

Y 

product  initially  defined  by 


(y,y)  *  (y.CpI  +  A)y)  ,  y,  y6«#(A), 

v  V  v 


being  selected  to  the  ^>1  +  A  is  strictly  positive. 
A  scalar  observation 


z(t) 


-  (y(t) ,  p  )  +  (X(t),  f  )  ,  p  S  Y,  f «  X, 


oo 


oo 


£  -  27  £  0  *  9“  27  9  0 

k*l  k  k  '  k-1  'k 


oo  .  2 


oo 


s lei <00.  £*.X)/»|< 

k«l  k  k*l  '  k  •  k 


oo, 


is  a  continuous  function  of  t  which  may  be  expressed,  assuming  it  to  be  real 
valued,  in  the  form 


cxj>  i«  t  oo  _  -iw.t. 
z(t)  *  2  e  K  «•  z7  ze  A 

k*l  k  k*l  k 


If  the  operator  A  is  perturbed  to  the  operator  A  +  S  A,  having  the  same 
properties,  but  with  eigenvalues  \  ,  "frequencies'*  c 0  +  ,  and  (non- 


normalized)  eigenvectors  0  +  $0  ,  ?0  6101  ,  and  if  the  same  observa- 

k  k  k  k 


tion  process  is  used,  there  will  result  the  modified  scalar  observation 
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00  i(u  *Sio)t  oo  -i((o  tSaJt 
<Ct)  -  ZT<e  k  k  ♦  £]?e  k  k 
’— ■*  •-  k«l  k 


k-1  k 


Through  the  use  of  discrete  Fourier  techniques  the  frequencies  CJ^+  (real 

exponents  if  X^<  0)  can  be  estimated,  though  with  some  additional  difficulty  if 

there  actually  are  X  <  0.  Then  the  identification  question  which  we  wish  to  ad- 

k  \ 

dress  is  this:  assuming  the  A  (hence  the  to  )  and  the  0  known  (i.e.,  the 

k  k  k  c 

operated  A  is  assumed  to  be  known  completely),  and  assuming  the  Co  +d  to  to  be 

C\  k  k 

"adequately  estimated",  which  means  that  the  dA  have  been  estimated,  can  we  det- 

k  c 

ermine,  at  least  in  first  order  approximation,  the  operator  perturbation  b  A? 

Before  going  on,  we  should  remark  at  this  point  that  any  identification  tech¬ 
nique  for  Sa  based  only  on  observations  ^(t)  must,  in  effect,  address  the  same 

problem,  for  the  i?  ,  S  are,  except  for  rates  of  decay,  etc.,  arbitrary  sequences 
yk  7k 

if  the  solution  of  the  system  is  unrestricted  and,  therefore,  the  only  real  infor¬ 
mation  available  consists  of  the  S/O  ,  equivalently  the  £X  ,  albeit  perhaps  in 

k  k 

an  indirect  form. 

The  first  order  analysis  (see  17  3,  1 13  I  for  more  details)  proceeds  as 

follows.  Accepting  the  equations 

(-(X  +£X  )I+(A+$A))(0  +  S0  )  =  0,  k  =  1,2,3,..., 
k  k  k  k 

we  take  the  inner  product  of  each  left  hand  side  with  the  corresponding  0  to 

c  ±  k 

obtain ,  taking  into  account  d  0  6.  t  0  3  , 

k  k 

(0  ,  SA0  )  *  $X  ,  k  »  1,2,3, ... 
k  k  X  k 

In  some  cases  the  equation  is  originally  given  in  the  form 
My  +  Ay  =  0 

with  M  bounded,  self-adjoint,  and  strictly  positive.  To  study  variations  M  ■> 

M  +  Sm,  A  remaining  fixed,  one  now  takes  the  0  to  be  the  solutions  of  the 

modified  eigenvalue  problem  (-Xm  +  A) (7  *»  0.  The  corresponding  equation  is  then 

k 


(0  ,  M0  )  SX  ♦  X  (0  ,  &M0  )  =  0,  k  =  1,2,3,... 
kkXkkk  kX 

Assuming  the  0  normalized  relative  to  M,  the  above  becomes 
k 

(0  ,  M0  )  *  -  £A  /  X  =  -  iogX  | 

k  k  k  k  » 

k 

which  we  can  use  in  the  same  way  to  infer  properties  of  the  perturbation  Sm. 
This  second  problem  is  usually  simpler  than  the  first  and  will  be  used  for  most 
of  the  examples  which  we  cite  here. 

Consider  the  more  specific  case  wherein  y  +  Ay  *  0  is  the  abstract  version 

i 

of  the  partial  differentail  equation 

m(x)y  +  £y  =  0,  t  £  0,  0  4  x  £  L, 
tt 


where  m(x)  is  the  mass  per  unit  length  of  the  medium  occupying  the  spatial  inter- 
val  10, L2  and  ,  the  elasticity  operator,  is  paired  with  boundary  conditions 


making  it  a  self-adjoint  operator.  If  the  normalized  eigenfunctions  are  0  (x), 

k 

k  *  1,2,3,...,  i.e., 

m(x)0  (x)0  (x)dx  »  a  > 

0  k  j  kj 


/' 

dO 


then  the  problem  of  determining  to  first  order  a  small  mass  variation  m(x)  leads 


to  the  infinite  set  of  equations 
•L  2 

®  (x)  0  (x)  dx  =  ■  $  /  /  ^  ,  k  *  1,2,3,..., 
Ok  k  k 


>V 

s: 


constituting  a  moment  problem  on  0  ^  x  4  L.  We  will  assume  that 

lim  $A/A  =  A 
k*oo  k  k 


exists.  Then,  since 

*L  2 

m(x)(0  (x))  dx  =  1,  k  =  1,2,3,..., 
0  k 


>  ***» 
/: 


replacing  £m(x)  by  Sm(x)  -  Am(x)  we  arrive  at  a  similar  problem  with  A  replaced 
by  zero.  Thus,  without  loss  of  generality,  we  may  assume 
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lim  £A  /  A  =0. 

k*oo  k  k 

The  rate  at  which  t\  /  \  tends  to  zero  will  ordinarily  determine  our  choice  of 

k  k  c 

function  space  for  m(x)  to  lie  in.  It  should  be  noted  that  if  we  think  of  c)m(x) 


as  a  mass  rotating  about  the  x-axis  at  a  distance  0  (x) ,  then  the  expression 
>L  2  k 

Sm(x)(0  (x))  dx 
0  k 


ALJ-UJ 

r. 


is,  in  fact,  the  moment  of  inertia  of  o m(x)  about  the  x-axis.  Thus  the  moment 

problem  in  fact  consists  of  finding  £m(x)  with  moments  of  inertia  -  §\  /2  A  , 

k  k 

k  *  1,2,3,...  (assuming  no  A  *  0,  of  course). 

k 

In  this  paper  we  will  confine  attention  to  the  two  questions  of  existence 
and  uniqueness.  For  the  existence  question  we  suppose  the  arbitrary,  sub¬ 
ject  to  having  an  approDriate  rate  of  decay  as  k*oo,  and  ask  if  a  solut- 

k  k 

ion  m(x)  exists  in  a  corresponding  function  space.  The  uniqueness  question  is 

posed  in  an  equally  self-evident  way.  It  turns  out  that  uniqueness  essentially 

never  obtains  without  further  restriction  on  the  form  of  Sm(x).  The  resolution 

of  the  existence  question  varies  from  one  example  to  another.  It  is,  in  reality, 

2 

the  question  of  the  linear  independence  of  the  functions  0  (x)  in  the  appropri- 

2  k  2 

ate  space,  which  we  usually  take  to  be  L  £o,L}.  If  the  functions  0  (x)  ,  k  -  1, 

k 

2,3,...,  are  strongly  independent,  i.e.,  there  exist  biorthogonal  functions  /l/(x) , 

2  & 
k  »  1,2,3,...,  also  in  L  £o,lJ,  such  that 

J-»L  2 

(x)  0  (x)  dx  =  *■* 

0  »  k  j 

then  we  have  a  formal  solution 


hat  _ 
kj  [0, 


k  *  j, 
k  j  j. 


m 


00  *  #Vxk)fkw 


whose  convergence,  depending  on  the  properties  of  the  y  (x) ,  can  be  investigated 

'  k 

subsequently. 
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There  is  .othing  about  the  problem  in  the  abstract  which  can  guarantee 
that  the  linear  independence  referred  to  above  obtains  in  general.  If  we  take 
m(x)  *  1,  L  *  2'ft'  ,  and 

)Cy  *  -y" 

with  periodic  boundary  conditions 

y(0)  *  y(2<r) 

'  2 

then  the  eigenvalues  are  \  =  0,  \  *  X  =  k  ,  k  *  1,2,3,...; 

o  2k  2k-l 

all  but  X  having  multiplicity  2  and  the  corresponding  eigenfunctions  are 
o 

0  (x)  *  cos (kx)  ,  -0  (x)  »  sin(kx)  ,  k  *  1,2,3,..., 

2k  2k- 1  VF 

while  the  simple  eigenvalue  0  is  associated  with  the  eigenfunction  0  fx)  -  1  . 

o  i 7Rr 

The  squared  eigenfunctions  are  then 

2  2 

1/xfr'  .  (l/1f)cos  (kx),  (l/1T)sin' (kx)  ;  k- »  1,2,3,..., 

2  2 

and,  since  (l/lrOcos  (kx)  *  (1/10  -  (l/IOsin  (kx),  this  set  of  functions  is 
certainly  not  linearly  independent  in  any  sense.  In  this  case  a  necessary  con¬ 
dition  for  solvability  of  the  moment  problem  is  clearly 

SA  *  o,  J  A  *-£A 

o  2k  2k-l 

and  a  general  existence  result  for  solutions  of  the  moment  problem  cannot  be 
obtained  for  arbitrary  frequency  perturbations  via  a  first  order  analysis. 

What  may  fairly  be  regarded  as  a  "prototype"  case  arises  again  for  m(x) 

»  1,  L  3  ,  and 

y  ■  »y i • 

but  with  separated  boundary  conditions 

a  y(0)  +  b  y*(0)  *  0,  a  yflT)  +  b  y'Otf)  =  0, 
o  o  1  1 

2  2  2  2 
a  ♦  b  j*  0,  a  +  b  jf  0. 
oc  11 


In  all  of  these  cases  the  eigenvalues  are  non-negative  and  simple,  the  zero 

eigenvalue  X  =  0  occurring  only  in  the  case  a  »  a  *  0.  There  are  four 
o  o  1 

principal  cases  to  consider: 

2 

Case  1:  a  *  a  =0.  Here  X  *0,  X  »k  ,k»  1,2,3,...,  and 

o  1  ok 

0  (x)  =  1//T  ,  0  (x)  =  (1/^21?)  cos  (kx) . 
o  k 

2  2 

Case  2:  b  /  0,b  /  0,  a  +a  /  0.  Here  the  eigenvalues,  all  posi- 

o  1  o  1 

tive,  have  the  form 

\  -k  ♦<?(!),  k  =  0,1,2,..,, 
k 

and  with 

CO  ■K  ->  +&  (l/k)  • 

k  k 

the  eigenfunctions  take  the  form 

0  (x)  *  y  cos(0>x  +  %  ), 
k  k  k  k 

where 

j)  *  u/i/2tb  +c?(i/w,  &  =  c9" ci/k) . 

k  k 

Case  3:  b  *  0,  b  f  0  (the  "mirror  image"  case  b  /  0,  b  =0  will  not  be 
o  1  ol 

treated  separately) .  Here  the  \  are  again  all  positive  and,  asymptotically, 

k 

X  -  (M0  +$(1), 

k 


while  the  eigenfunctions  and  frequencies  are  given  by 

0  (x)  »  -j/sinCG}  x),  k  *  1,2,3,..., 
k  k  k 

Cl)  =  (k-y  +  0 (l/k) ,  i)  =  (1/-J2V)  +  #(l/k),  k  =  1,2,3,... 
k  k 

Case  4:  b  »  b  »  0.  Here  \  *  k  ,  k  *  1,2,3,...,  and 
o  1  k 

0^(x)  «  (1/1/21? ) sin (kx) ,  k  «  1,2,3,...  . 
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*■ 


r 


As  noted  earlier,  in  all  cases  we  may  assume  the  variation  Sm(x)  replaced 
by  m  +  £m(x),  where  m  is  a  constant  and 

£m(x)dx  =  0. 

'0 

2  2 

Since  (l/2fl')cos  (toe)  =  (1/21?)  -  (l/21t/)sin  (kx)  and  £m(x)  has  been  assumed 

orthogonal  to  constants.  Cases  1  and  4  may  be  treated  together,  without  loss  of 

generality  as  Case  1.  Here 

2 

(1/210 cos  (kx)  =  (l/4ff)(l  +  cos (2kx) ) , 
and  again  making  use  of  the  above  orthogonality  property,  the  moment  problem 
reduces  to 


cos(2kx)dx  *  -4TT  S  X  /k  ,k*  1,2,3,...  , 

k 


the  case  corresponding  to  X  “  0  already  having  been  taken  care  of,  assuming 

SA  =  0.  The  functions  cos(2kx)  are,  of  course,  orthogonal  on  [o ,1Tjand, 
o 

hence,  independent.  Adjoining  1,  sin2kx,  k  »  1,2,3,...,  and  normalizing  we  have 

2 

an  orthonormal  basis  for  L  £o,ft'J.  The  unique  solution  of  the  moment  problem 
orthogonal  to  constants  and  even  in  the  sense 
$m(tr  -  x)  «  Sm(x) 

is,  formally, 

oo  2 

m(x)  »  -2 T  27(^X  A  )cos(2kx)  , 
k«l  k 

2  -  c\  (3/2)-  S) 

and  is  convergent  in  L  {.0 ,  ffj  as  long  as  a  A  <*  C/  (k  ). 

k 


If  the  even-ness  requirement  is  removed,  the  solution  is  highly  non-unique,  since 
one  can  add  to  Sm(x)  any  function  m(x)  in  L  [o, If ]  which  is  odd  in  the  sense 
ra(TT -  x)  =  -m(x) . 

Case  3  is  mathematically  nontrivial  but,  fortunately,  the  work  has  already 
been  done  for  us .  It  is  known  from  the  work  of  Paley  and  Wiener  1 12  I ,  Levinson 
I  9  I,  Schwartz  I  15],  and  many  others,  that  the  functions 
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£in(2o>x),  cos (2 u>  x) ,  ■  (k-y  +  $(l/k),  k  =  1,2,3,... 

k  k  k 

2 

form  a  Riesz  basis  for  L  ,  provided  the  CO  are  distinct,  as  in  the 

k 

case  here.  Among  other  things  this  implies  the  existence  of  a  unique,  biortho- 

gonal ,  dual  basis  p  (x),  q  (x),  with 

k  k  ^ 

Xsin(2ttfx)p  (x)dx  =  /  cos(2tc>x)q  (x)dx  =  $ 

*  k  j  ^0  k  j  kj, 

/  sin(2a>x)q  (x)dx  *  /  cos(2£Ox)p  (x)dx  =  0. 

«/  0  k  j  wn  k  j 

2 

Each  function  f  in  L  [b^Jhas  unique  expansions 

oo  '  oo  A  A 

fCx)  *  «Zj(c  -sin (2 CM  x)  +  d  cosC2u>  x))  =  ^(c  p  (x)  +  d  q  (x)), 
k»l  k  k  k  k  k*l  k  k  k  k 

2 

convergent  in  L  £o,HfJ  with 


c  *  /  p  (x)f(x)dx,  d  .  ■  f  q  (x)f(x)dx, 
k  v0  k  k  "0  k  ^ 

a  r*  a  r* 

c  *  /  sm(24c;  x)f(x)dx,  d  -  J  cos(2 CO  x)f(x)dx 
k  "0  k  k  "0  k 


and,  for  certain  positive  numbers  R,  r 


r  Ilf  ||  2  ^  J§(  (c|  +  U  I  )  ^  Rff  I  2 

L  fo/tr'J  k=l  k  k  L  [O,*! 

(l/R)  Ilf  II  2  *  X(  |c  l  +  Id  |  )  ^(1/r)  ||f  J  2 

L  (O/tf]  k=l  k  k  L 


The  moment  problem 
/’TT 


f  $ m(x)0  (x)  dx  «  -  S  X  /  A  »  k  *  1,2,3,... 
c/O  k  k  k 


in  this  case  becomes  a, 

2  /*"  *  2  2 
V  /  £.u(x)sin  (.  Ct7  x)dx  *  -  S'X  /Ck  +  0(1)),  k  =  1,2,3,...  , 
k  vr  k  k 

and,  since  Sm(x) ^is  assumed  orthogonal  to  constants,  this  is  the  same  as 
f  Sm(x)cos(2Cc>  x)dx  »  2$  X  /(-^  (k  +  (?(!))) 

<J(l  lr  V  V 


The  unique  solution  lying  in  the  space  spanned  by  the  functions  q  ,  j  =  1, 

3 

2,3,...  is 

oo  .  2  2 

m(x)  =  2  S(SA/(yCj  +  &  Ci))))q  Cx) 

j  *  1  j  j  J 

convergent  just  in  case  the  sequence  of  coefficients  is  square  summable.  Again 

this  solution  is  non-unique;  one  may  add  any  convergent  linear  combination  of 

the  p  (x)  and  still  have  a  solution  of  the  moment  problem.  Here  the  displayed 

solution  cannot  be  characterized  as  odd  or  even  relative  to  the  center  point  If  /2 

of  the  x-interval  because  the  functions  cos(2  6Jx),  sin(2*y  x)  have  no  such  prop- 

k  k 

erties,  except  in  some  asymptotic  sense  as  k  »  oo. 

Finally,  Case  3  is  the  most  complex  because  here 
2  2  2  2 
0  (x)  *7/  cos  (fit)  x  +  S  )  *  l)  C(1  +  cos (2 CO  x  +  2  S  ))/2)  . 
k  k  k  k  k  k  k 

The  constant  term  can  be  discarded  for  the  same  reasons  as  before  and,  since 

cos(2co  x  +  2  J  )  =  cos(2  §“  )cos(2 (o  x)  ••  sin(2  £  )sin(2£<>  x), 
k  k  k  k  k  k 

the  relevant  moment  problem  becomes  ^ 

r«  (% 

J  dm(x)cos(2Cc?  x)dx  +  p  J  d  m(x)sin(2 Ct)  x)dx 

«  -(2  )/(V  ((k-32)2  +  1))), 

k  k 

where,  as  k  4  00, 

2 

ol  »  cos (2  £  )  =  I  +  0  (1/k  ), 
k  k 

&  «  -sin  (2  S  )  =  <^Cl/k). 

Ik  k 


Perhaps  the  most  natural  solution  here  is 


Sn(x) 


JJcSX/Cz/CCj-J*)2  +^d))))foCq  (x)  ♦  Ap  Cx)] 
j«l  j  j  j  j  'j  j 

with  convergent  linear  combinations  of  the  functions  5  q  (x)  -  p  (x) 

\j  j  j  j 


constituting  the  space  of  functions  which  may  be  arbitrarily  added  to  £m(x),  as 
displayed,  without  changing  its  property  of  being  a  solution  of  the  moment  problem. 

If  we  maintain  m(x)  =  1  and  suppose  that  the  equation  m(x)y  +^y  *  0  is,  in 

tt 

reality, 

y  -  (p(x)y  )  =  0, 
tt  x  x 

then  perturbations  p00  =  1  +  £ p (x)  about  the  nominal  p  (x)  =  1  lead,  as  noted 

o 

previously,  to  the  equations 

Vir  2 

/  0  (x)(£p(x)0'(x))  *dx  =  /  apCx)C0'Cx))  dx  =  d\  ,  k  «  1,2,3,..., 

J o  k  •  k  "0  k  k 

which  may  be  analyzed  in  much  the  same  way  as  the  problem  for  ?m(x)  treated 
in  detail  above  with  comparable  solutions  and  non-uniqueness  aspects  as  were  de¬ 
veloped  there . 

To  conclude  this  section  on  identification  from  frequency  data,  let  us  con¬ 
sider  the  case  of  the  Euler-Bemoulli  beam,  which  corresponds  to 
m(x)y  +  X  y  a  0 


Xy  *  (EUx)y  ) 

xx  xx 

and  appropriate  boundary  conditions,  for  example,  the  cantilever  beam  (the  clamped/ 
free  case)  corresponds  to 

y(0,t)  =  y  (0,t)  *»  0, 
x 

y  (L,t)  =•  (El(x)y  )  *  0. 

rx  xx  x  x  *  L 

If  we  begin  with  m  (x)  3  i,  EI(x)  =  1,  L  »  ,  and  consider  m(x)  =  1  +  $m(x) 

o 

r* 

for  small  om(x)  we  are  again  led  to  the  moment  problem 

/  $m(x)(0  (x))  dx  *  -  ^  /  \  . 
va  k  k  k 


f 


As  we  have  seen  earlier,  if  in  this  case  we  set  CO  =  ■)  \  ,  the  CO  are  the  posi- 

k  k  k 

tive  roots  of 

h  h 

1  +  cos  (It?  'R')cosh(&;  /it)  *  C 


h  h  -1 

cos  (a?  It')  =  -CcoshC^  Tt )) 


leading  to 


h  .  -kit* 

6/  =  (k-%)  +  0  (e  )  . 
k 


Then  from  the  equation 

h  h  h  h 

c(-sin(OJ  It*)  -  sinh(ii>  1h))  +  d(-co s(CJ  it)  -  cosh(&;  It))  -  0 


it  is  clear  that  for  (0  -  0)  we  must  have 

k 


c  *  c  a  tO  +  y  e  ^ 
k  k  k  , 

d  *•  u  »  -j/  +  §  e 
k  k  k 


where,  if  the  eigenfunctions  0  (x)  are  normalized, 

k 


lim  x)  * 
k*oo  k 


for  some  positive  constant  %/  and  the  Y' ,  §  are  uniformly  bounded.  Thus 

.  k  k  ' 

0  (x)  =  (V  +  y  e  k  )(sin(&Ac)  +  e  k  ) 
k  k  k  k  ,  . 

^  -o?ir  h  -co\  -^t('r-x) 

+  %  e  )(cos(«  x)  -  e  k  )  +  (Y  -  §“  )e 
k  k  k  k  k 


From  this  it  may  be  seen  that  0  (x)  is  a  linear  combination,  with  bounded  coef- 

k 

ficients,  of  the  functions  ,  . 

Jr  '  -^gx 

sin(w  x),  cos(cjsx),  e  ,  e 
k  k 


In  this  case  we  know  from  the  theory  of  the  operator  itself  that  the  0  (x  ) 

2_  -  k 
are  orthonormal  in  L  [0,1*4 ,  and  hence  independent,  but  no  theorems  are  readily 


ru  wdj.w^jLU  ■  .ilw,iuuhi  wi  tf -^gw-wmB*  gjhVJUPimn? 


4\ 


available  concerning  the  independence  properties  of  the  functions  (0  (x))  , 

k 

which  constitute  the  defining  kernels  for  the  moment  problem  equivalent  to  the 


first  order  frequency -based  identification  problem  for  $n(x).  Interest  in 
this  problem  will  thus  lead  to  further  study  of  tj re  independence  properties 
of  complex  exponentials  and  related  functions  (see  1 14  J ,  [  16  ]  ■) . 

The  results  noted  above  for  the  vibrating  string  all  have  the  characteristic 
that,  in  terms  of  density  of  the  exponents  of  the  exponential  functions  involved, 
only  half  as  many  equations  are  available,  based  on  observed  frequency  perturbat¬ 
ions,  as  would  be  required  to  determine  £m(x),  or  'bp(x) ,  completely.  If  a 
boundary  control  force  can  be  implemented,  consistent  with  the  boundary  condition 


y  (tf ,t)  »  u(t), 
x 

for  example,  we  may  obtain  the  complementary  set  of  equations  by  determining  u 


via  boundary  feedback .  For  example,  with 


u(t)  *  -Y y  (If  ,t) , 

t 


another  set  of  perturbed  spectra,  and  related  moment  equations,  can  be  obtained, 
corresponding  to  the  boundary  condition 


y  (fT ,t)  +  y  (fl'.t)  *  0 
x  t 


and  these,  along  with  those  obtained  for  u(t)  *  0,  can  be  shown  to  be  sufficient 
to  determine  $m(x),  or  %  p (x) ,  completely.  This  is  in  agreement  with  the  more 
global,  but  theoretical,  results  c-f  Borg  11  1,  18  I. 

It  is  likely  that  j.n  most  applications  the  identification  problem  will  be 
less  complicated  than  the  one  we  have  treated  here.  In  many  cases  we  may  exp<j*.t 
ohat  it  will  bi  known  in  advance,  for  the  mass  perturbations  £m(x),  for  example, 
that 


m(x)  *  Ln  (x)  +  um  (x)  +  . ..  ♦  urn  (x). 


Vi 


•  urn  l.xj  +  ...  +  um 
1  /  2  2  /  r  r 


3$ 


where  m  (x) ,  m  (x) , . . .  ,m  (x)  are  functions  known  in  advance .  This  would  be  the 
12  r 

case,  for  example,  if  the  individual  terms  represented  fuel  or  payload  of  un¬ 
known  mass  but  stored  in  known  configurations  on  the  elastic  body  in  question. 

In  such  cases  the  moment  equations  in  fact  overdetermine  the  coefficients  jji 
to  be  identified;  some  sort  of  least  squares  solution  has  to  be  implemented, 
which  leads  to  an  unambiguous  result  provided  that  no  nontrivial  linear  combi- 

2 

nation  of  the  m  (x)  is  orthogonal  to  all  of  the  squared  eigenfunctions  ($  (x))  . 

k  k 
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